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We update our previous results by increasing the luminosity, the efficiency, and for the DJ" the number of tags. 
We determine f D+ = (205.8 ± 8.5 ± 2.5) MeV, and an interim preliminary value of / + = (267.9 ± 8.2 ± 3.9) 
MeV, where both results are radiatively corrected. We agree with the recent most precise unquenched Latticc- 
QCD calculation for the D + , but are in disagreement for the D£ . Several consequences are discussed, including 
the possibility of physics beyond the Standard Model. 
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1. Introduction 

Purely leptonic decays of heavy mesons proceed in 
the Standard Model (SM) via a W + annihilation dia- 
gram shown specifically for D + — > in Fig. 1. The 
strong interaction effects are parameterized in terms 
of the "decay constant" for the D + meson f D + . The 
decay width is given by 
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Figure 1: The decay diagram for D + — > l + v. 
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where Gf is the Fermi constant, M D + is the D + mass, 
mi the final state charged-lepton mass, and V c d is a 
CKM matrix element, taken equal to V us . Thus, in 
the SM measurement of purely leptonic decays allow 
a determination the decay constant, f D + of the D + 
meson and similarly fjj a of the meson. 

Meson decay constants in the B system are used to 
translate measurements of BB mixing to CKM matrix 
elements. Currently, it is not possible to determine 
Jb accurately from leptonic B decays, so theoretical 
calculations of Jb must be used. Since the B® meson 
does not have £ + v decays, it will never be possible to 
determine fs B experimentally, so again theory must 
be relied upon. If calculations disagree on D mesons, 
they may be questionable on B mesons. If, on the 
other hand new physics is present, it is imperative to 
understand how it effects SM based predictions of the 
B decay constants. Decay constants can be calculated 
using Lattice-QCD techniques. Recently, Follana et 
al. using an unquenched lattice technique predicted 
f D + = (207 ± 4) MeV and f Ds = (241 ± 3) MeV. [1] 



In these analyses we exploit the reactions e + e — » 
D-D+, and e+e" -> D*~D+ or DjD* s +. The D+ 
is studied at 3770 MeV using 818 pb _1 . D+ is 
studied at 4170 MeV, using 400 pb _1 for the 
and t + v, t + — > 7T+T7 final states, and 300 pb _1 for 
t + — > e+vv. (Eventually CLEO will present results 
using 600 pb _1 .) 



2. D+ -»■ l+v 

We use a "double tag" technique where one D ± is 
fully reconstructed and the oppositely charged D can 
then be found even if there is a missing neutrino in the 
final state [2]. For notational convenance, the D~ is 
referred to for the fully reconstructed tag, although 
D + states are also used. To reconstruct D~ tags 
we require that the tag candidates have a measured 
energy consistent with the beam energy, and have a 
"beam constrained mass," msci consistent with the 
D~ mass, where m BC = v^bcam = EiPiFi E hcam 
is the beam energy and i runs over all the final state 
particles three- momenta. Fig. 2 shows the tubc dis- 
tribution summed over all the decay modes we use for 
tagging. Selecting events in the mass peak we count 
460,055±787 signal events over a background of 89,472 
events. 

To search for signal events we look for events with 
one additional track with opposite sign of charge to 
the tag. The track must have an angle >25.8° with 
the beam line. We separate these events into two cat- 
egories. Case (i): those which deposit < 300 MeV 
of energy in the calorimeter, characteristic of 99% of 
muons, and case (ii) those which deposit > 300 MeV, 
characteristic of 45% of the pions, those that happen 
to interact in the calorimeter. 

We look for D + — > n + v by computing the square of 
the missing mass 
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where Pn- is the three- momentum of the fully recon- 
structed D~ , and B„+(p„+) is the energy (momen- 
tum) of the [i + candidate. The signal peaks at zero 
for fi + is and is smeared toward more positive values 

for T + V, T + — > TT + V. 
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Figure 2: The beam-constrained mass distributions 
summed over D~ decay candidates in the final states: 

K + TT~TV~ , K + TV~TT~-K° , KsTV~ , Ks TV ~ 7T ~ 7T + , KsTV~TV° and 



The fit to the case (i) MM 2 distribution is shown 
in Fig. 3 contains separate shapes for signal, 7r + 7r°, 

K 7T + , T + l> (t + — 



the 



tt + D), and a background shape 
describing three-body decays. Here we assume 
SM ratio of 2.65 for the ratio of the "- + 
component and constrain the area ratio of these 
components to the product of 2.65 with B(t + — > 
7r+^)=(10.90±0.07)% [5] and the 55% probability that 
the pion deposits <300 MeV in the calorimeter. We 
veto events with an extra neutral energy cluster > 250 
MeV. This removes most 7r + 7r° events; the residual 
background is fixed at 9.2 events. The normalizations 
of the signal, K tt + , and 3-body background shapes 
are allowed to float. 

The fit yields 149.7±12.0 n + v signal events and 25.8 
t + v, t + — > tt + D events (for the entire MM 2 range). 
We also perform the fit allowing the t + v, t + — > ir + D 
component to float. Then we find 153.9±13.5 fi + v 
events and 13.5±15.3 t + v, t + — > 7r + P events, com- 
pared with the 25.8 we expect in the SM. Performing 
the fit in this manner gives a result that is indepen- 
dent of the SM expectation of the D + — > t + v rate. To 
extract a branching fraction, in either case, we sub- 
tract off 2.4±1.0 events found from simulations and 
other studies to be additional backgrounds, not taken 
into account by the fit. 

Wefind£(L>+ -> n + v) = (3.82±0.32±0.09) x 1(T 4 . 
The decay constant f D + is then obtained from Eq. (1) 
using 1040±7 fs as the D+ lifetime [5] and 0.2256 as 
\Vcd\- Our final result, including radiative corrections 
is 




f D + = (205.8 ± 8.5 ± 2.5) MeV . 



(3) 
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Figure 3: Fit to the MM 2 for case (i). The points with er- 
ror bars show the data. The black (dashed) curve centered 
at zero shows the signal fi + u events. The dot-dashed (red) 
curve that peaks around 0.05 GeV 2 shows the D + — > t + v, 
t + — v tx + u component. The solid (blue) Gaussian shaped 
curve centered on the pion-mass squared shows the resid- 
ual tt + tv° component. The dashed (purple) curve that falls 
to zero around 0.03 GeV 2 is the sum of all the other back- 
ground components, except the K 7r + tail which is shown 
by the long-dashed (green) curve that peaks up at 0.25 
GeV 2 . The solid (black) curve is the sum of all the com- 
ponents. 

3. D+ £+l> 

In the D s case we have to take into account the 
additional photon from D* — > jD s , since we use 
e + e~ — > D*D S events. We first examine the invariant 
masses (see Fig. 4). Then to select the appropriate 
D~ tag and 7 candidate, we compute the square of 
the missing mass opposite the selected tag and can- 
didate 7's, which peaks at the Df mass for correct 
combinations, where 



MM* 2 = (£ CM - E Ds - E n f - (p CM - p Ds - P 7 ) 2 , 

(4) 

here Eqm (pcm) is the center-of-mass energy (mo- 
mentum), Ed, (pd b ) is the energy (momentum) of the 
fully reconstructed D~ tag, and E 1 (p 7 ) is the energy 
(momentum) of the additional photon. In performing 
this calculation we use a kinematic fit that constrains 
the decay products of the D~ to the known D s mass 
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and conserves overall momentum and energy. 
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Figure 4: The invariant mass distributions summed over 
D~ decay candidates in the final states: K + K~ n~ , 
K$K~ , nn~; n — » 77, n'n~; n 1 — > n + n~w, n — > 77, 
4>P~\ 4> -> K+K-, p- -> ir-TT°, iy + tt~iy-, K*-K*°; 



and n'n~; n' 
background. 



K*° — > .K" + 7r , ?7p ; 77 — > 77, p — * ty tt u , 
* 7r + 7r~7. The curves represent signal and 



The MM* 2 distributions from the selected D~ event 
sample are shown in Fig. 5. We fit these distributions 
to determine the number of tag events. 
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Figure 5: The MM* 2 from a 7 plus tag candidate. The 
curves represent signal and background components. 

After selecting the tags we then find events with a 
single oppositely charged track to the tag and compute 



MM = (E C m - /'./, - E 
— (pcm — Pd s 



1 

Pi 



(5) 



We make use of a set of kinematical constraints and fit 
each event to two hypotheses one of which is that the 
tag is the daughter of a D*~ and the other that 
the D* + decays into 7.D+, with the Df subsequently 
decaying into The MM 2 distributions from data 



are shown in Fig. 6 where we have summed cases (i) 
and (ii). After fixing the ratio of t + v/ ji + v to the SM 
value we find f D + = (268.2 ± 9.6 ± 4.4) MeV. 
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Figure 6: The MM 2 distribution. The dashed (grey) 
Gaussian shaped curve peaked near zero is the p + v com- 
ponent, while the dashed (purple) curve that rises sharply 
from zero and then flattens out shows the t + u component. 
The two lines are background components. The solid curve 
shows the sum. 



We can also use the decay mode r + — > e + vv. This 
result has already been published. [6] The technique 
here is to use only three tagging modes: <pTr~ , K~K*° 
and KgK~, to ensure that the tags are extremely 
clean. Then events with an identified e + and no 
other charged tracks are selected. Any energy not 
associated with the tag decay products is tabulated. 
Those events with small extra energy below 400 MeV 
are mostly pure — > t + v events. After correct- 
ing for efficiencies and residual backgrounds we find 



f D + = (273 ± 16 ± 8) MeV. 



4. Conclusions 

The preliminary CLEO average is f D + = (267.9 ± 
8.2 ± 3.9) MeV (radiatively corrected). Averaging 
in the Belle result [7] f D + = (269.6 ± 8.3) MeV, 
which differs from the Follana et al. calculation [1] 
by 3.2 standard deviations, while the result for fjj+ = 
(205.8±8.5±2.5) MeV is in good agreement. This dis- 
crepancy can be explained either by New Physics [4] or 
casts suspicion on the theoretical prediction. As sim- 
ilar calculations are used for /b s //b, we need worry 
about them, or the effects of New Physics on this ra- 
tio. 
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